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Abstract 

The vertex separator problem (VSP) in an undirected connected graph G — 
{V, E) asks for a partition of V into nonempty subsets A, B, C such that |C| is 
minimized subject to there is no edge between A and B, and maa;{jA|, \B\} < 
[3{n) with P{n) integer such as 1 < /3(n) < n. We investigate the (VSP) in a 
polyhedral way, starting from the formulation given by E. Balas and C. De Souza 
m . We introduce new efficient valid inequalities. Our algorithm is based on some 
of the inequalities discussed here. The power of thèse inequalities is illustrated by 
the computational tests, efficiently improving the Balas and De Souza's results t2J . 

Résumé 

Le problème du séparateur (VSP) dans un graphe connexe non orienté G — 
{V, E) est celui de la détermination d'une partition de V en trois classes non vides 
A, B, G telle que jC| soit minimum sous les contraintes qu'il n'existe aucune 
arête entre A et B et que maa::{|yl|, \B\} < /3(n), où /3(n) est un entier tel que 
1 < /3(n) < n. Le problème est étudié d'un point de vue polyédrique à partir 
de la modélisation de E. Balas et C. De Souza 1 1 1 que nous avons enrichie. En 
particulier, nous introduisons de nouvelles inégalités valides qui se révèlent d'une 
grande efficacité dans la résolution de certaines instances difficiles du problème. 
Notre algorithme est basé sur certaines des inégalitées présentées ici. Les résultats 
obtenus améliorent très sensiblement ceux de Balas et De Souza [2- 
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1 Introduction 

Dans un graphe connexe non orienté, un séparateur est un sous-ensemble de som- 
mets ou d'arêtes dont la suppression déconnecte le graphe. La recherche d'un sé- 
parateur de poids minimum équilibré (dont la suppression isole deux sous-ensembles 
de sommets de tailles équivalentes, au plus 2n/3) intervient dans de nombreux 
problèmes. 

Dans le domaine des réseaux et télécommunications, un tel séparateur est vu 
comme un goulot d'étranglement lorsque le réseau est représenté par un graphe. 
Lors de l'évaluation de la capacité du réseau et de sa robustesse, il est utilisé pour 
la recherche de bornes inférieures et l'identification des noeuds sensibles. 

Dans le domaine des algorithmes opérant sur les graphes, la connaissance des 
séparateurs de poids minimum équilibrés est fondamentale, en particulier pour 
ceux reposant sur le paradigme diviser et conquérir . 

Pour exemple, on peut citer les stratégies de résolution des systèmes linéaires 
creux de grande taille qui font appel à ce type d'algorithmes. La taille des prob- 
lèmes impose une puissance de calcul croissante et le recours aux calculateurs 
parallèles est un passage obligé. Ainsi, la factorisation des matrices creuses est 
particulièrement adaptée à la programmation parallèle. Plus la matrice est creuse, 
plus l'on est susceptible d'avoir de parallélisme. Il faut conserver le plus possi- 
ble l'aspect creux de la matrice durant la progression de la factorisation. Dans 
l'exemple de la factorisation de Cholesky, méthode directe de résolution des sys- 
tèmes carrés, symétriques définis positifs, la première étape de la factorisation 
consiste en une renumérotation des lignes de la matrice initiale pour obtenir une 
matrice de Cholesky la plus creuse possible |6|. La recherche de la meilleure per- 
mutation est un problème complexe nécessitant le recours à des heuristiques. Les 
algorithmes de renumérotation (degré minimal, dissection emboîtée) reposent sur 
les principes de la théorie des graphes. L'algorithme de dissection emboîtée repose 
sur le principe diviser et conquérir (divide and conquer). Il considère que la ma- 
trice initiale est irréductible (i.e. le système ne peut se décomposer en deux sous- 
systèmes indépendants) et qu'elle est donc associée à un graphe d'intersection con- 
nexeQ Le but de la dissection emboîtée est de faire interagir le moins possible des 
parties relativement indépendantes de la matrice. Pour celà, en partant du graphe 
initial G, on cherche un séparateur qui permettra la décomposition de G en deux 
sous-graphes qui pourront être traités en parallèle. La procédure est appliquée 
récursivement aux sous-graphes. Le parallélisme mis en évidence est tributaire de 
la taille des séparateurs : plus ils sont petits, meilleur est le parallélisme. Cet al- 
gorithme se décline selon plusieurs variantes (par exemple, la dissection emboîtée 
spectrale {Spectral Nested Dissection) (SND) Il II ) mais repose toujours sur la 
recherche de bons séparateurs. 

Dans le domaine du traitement automatique de la langue écrite, le concept de 
séparateur est utilisé en classification par Berry et al dans les grands graphes de 
termes creux f3\. Les sommets du graphe sont les termes du domaine consid- 
éré. Les arêtes du graphe représentent les relations syntaxiques entre les termes. 
L'objectif est d'obtenir, à partir de corpus textuels, des classes susceptibles de 

'Le graphe d'intersection d'une matrice A, noté G{A), a un sommet pour chaque colonne de A et une 
arête entre une paire de sommets si le produit scalaire des colonnes correspondantes est non nul. 

^syntagmes nominaux de plusieurs mots susceptibles de désigner, par leurs propriétés syntaxiques et 
grammaticales, un objet ou une notion du domaine. 
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représenter des thématiques. On procède par classification iiiérarchique 1121 mais 
certaines classes restent trop grosses pour être interprétables par l'utilisateur Elles 
sont alors décomposées en utilisant des séparateurs minimaux complets. La qualité 
de la décomposition est liée à la taille des séparateurs et à l'équilibre entre la taille 
des sous-classes obtenues. 

Dans le domaine de la bioinformatique, les séparateurs sont recherchés dans 
les graphes de grille modélisant la structure des protéines. Fu et al |7 1 explorent la 
résolution du problème des configurations d'énergie minimale qui doivent rassem- 
bler les sommets non consécutifs occupant des positions voisines. 

Formellement, le problème du séparateur (VSP) peut être énoncé de la façon 
suivante : 
Données : 

• G — {V, E) un graphe connexe non orienté avec \ V\ = n, 

• I3(n) entier tel que 1 < I3{n) < n, 

• Ci coût associé à chaque sommet i £ V. 

Objectif : 

Trouver une partition {A, B, C} de V telle que : 

E ne contient aucune arête avec i £ A ,j £ B, (1) 
max{|A|,|B|} < /3(n) , (2) 

Cj est minimum (3) 

jeo 



A C B 




Exemple de séparateur 



Le problème du séparateur (VSP) est un problème NP-difficile (4). 

Le plus ancien résultat concernant les séparateurs est obtenu en 1869 par Jor- 
dan. Il démontre alors que pour chaque arbre, un seul sommet le sépare en deux 
parties contenant chacune moins des deux tiers de ses sommets. Lipton et Tra- 
jan 1101 prouvent que tout graphe planaire possède un séparateur de taille au plus 
0{^/n) qui le sépare en deux parties contenant chacune moins des deux tiers de 
ses sommets. Dans le cas oîi G est planaire et /3(n) = ils établissent qu'un 
séparateur de taille bornée par 2\f2yfn peut être trouvé en un temps 0(n). Bui, 
Fukuyama et Jones 15J montrent que dans le cas des graphes planaires, la détermi- 
nation d'un séparateur optimal est un problème NP-difficile. 

En 2005, Egon Balas et Cid De Souza réalisent la première étude polyédrale du 
(VSP). Ils examinent le polytope des séparateurs, enveloppe convexe des vecteurs 
d'incidence des partitions \^A, B, C} de V vérifiant (î} et (0. Leurs recherches 
sont centrées sur la relation entre séparateurs et dominants. (Un ensemble dom- 
inant S d'un graphe G est un sous-ensemble de sommets de G tel que chaque 
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sommet de G est soit dans S, soit voisin d'un sommet de S). Ils donnent une for- 
mulation du problème à l'aide d'un programme linéaire mixte 1 1 1 et développent 
parallèlement un algorithme de branch-and-cut basé sur les inégalités issues des 
relations entre séparateurs et dominants (2). Les recherches de Balas et De Souza 
ont été notre point de départ. 

Notre travail repose également sur une approche polyédrale du (VSP). La for- 
mulation de Balas et de Souza utilise des variables sur les sommets du graphe G. 
Nous avons adjoint plusieurs inégalités valides à celles présentées. D'une part, 
nous avons considéré les chaînes entre chaque paire de sommets non ajacents du 
graphe G, ce qui nous a permis une première minoration du cardinal de tout sé- 
parateur associé à G. D'autre part, nous avons reformulé le problème en rajoutant 
des variables sur les arêtes. 

Dans la seconde partie de ce document, nous exposons l'approche polyédrale 
du problème. Nous étudions le cas particulier du VSP où a G A et 6 G B sont 
fixés. La dimension du polyèdre associé est donnée et des inégalités valides sont 
introduites. La troisième partie développe une modélisation du problème utilisant 
des variables sur les arêtes du graphe. Notre implémentation et les résultats que 
nous avons obtenus sont présentés dans la quatrième partie. La dernière partie 
expose nos conclusions et les perspectives de nos travaux. 
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2 Le polyèdre des séparateurs 

Pour un graphe G = {V, E) donné, on considère le cas particulier du VSP oîi deux 
sommets a et fe non adjacents sont donnés et l'on cherche une partition {A, B, C} 
satisfaisant (Tll et l|2j avec a £ A, b £ B eX telle que |C| soit minimum. Ce cas 
particulier sera nommé problème du ab-séparateur . Dans cette partie, on étudie le 
polyèdre associé au problème du ab-séparateur. En effet, la résolution du VSP se 
ramène à la résolution d'au plus sous-problèmes de ce type. 

Les sommets a et fo non adjacents étant fixés, le vecteur d'incidence d'une par- 
tition {A, B, C} vérifiant (î) et l[2j avec a G A elb £ S est donné par X G 

{0,l}2(n-2) : 

Xia = l<^i£ A, Xib = l<:^ie B,\lieV\{a,b} 

On désigne par Pab le polyèdre associé au problème du ab-séparateur, défini par : 
Pab = Conv{X G R^(""^' , 3 un ab-séparateur {A, B, C} dont X est un vecteur 
d'incidence }. 

Si X est un vecteur d'incidence d'un ab-séparateur {A, B, C}, alors X vérifie : 



Xia + Xjb<l yii,j)€E (4) 

Xja + Xib<l \f(i,j)eE (5) 

Xia + Xib<l yieV\ {a, b} (6) 

n — 2 

Y,Xia< 0{n)-l \fieV\{a,b} (7) 

1=1 

n-2 

^ < /3(n) - 1 Vi G V \ {a, fc} (8) 

1=1 

x,aÇi{0,l} ViGV\{a,fe} (9) 

x,b>0 yi€V\{a,b} (10) 



Les contraintes ([4J et (O sont valides du fait de l'absence d'arêtes entre les 
ensembles A et B. La contrainte ([6} est valide car un sommet ne peut appartenir 
à la fois a A et à B. Les contraintes ([71( et ([8} sont valides car les ab-séparateurs 
vérifient (O. 

Balas et De Souza (T| ont montré que 

Pab = Conv |x G r2("-2) , X vérifie @ à ([ÎÔ| j. 

2.1 Dimension du polyèdre Pab 

Un polyèdre P est de dimension k si le nombre maximum de points affinement 
indépendants de P est A: + 1. Si P Ç R'' alors dim{P) < k. 

Soient V{a) C et V{b) C V les ensembles de sommets adjacents respective- 
ment à a et à 6. 

Théorème 2.1 

dim(P^b) = 2(n - 2) - (|y(a)| + 
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Preuve : 

★ On obtient rapidement dim{Pab) < 2(n - 2) - (|l^(a)| + |1^(&)|). En effet, 
pour tout sommet i de V{a) (resp. de V{b)), on a Xib = 0, VX £ Pat (resp. 
Xia ~ 0, VX £ Pab)- La description minimale de Pat contient donc au moins 
|F(a)| + |1^(b)| égalités, d'oîi la majoration de dim{Pab)- 

L'égalité est obtenue par récurrence sur le nombre de sommets du graphe : 

★ Pour n = 3, on a Pat ~ {(0, 0)} d'oîi dim{Pab) ~ 0. 

Or 2(n - 2) - {\V{a)\ + \V{b)\) = 0, l'égalité est donc vérifiée. 

* Supposons l'égalité vérifiée pour tout graphe dont le nombre de sommets est 
inférieur ou égal à n — 1 et soit G — {V, E) tel que \ V\ = n, {a, b} C V. 

• Cas 1 : 3uo G V \ {a, b} tel que vo ^ V{a) U V{b). 

On suppose sans perte de généralité que vo porte le numéro {n — 2) dans G. 
Considérons G', le sous-graphe de G induit par V^' — V\{vq}. SoitP^j le polyè- 
dre des ab-séparateurs associés au graphe G'. 

D'après l'hypothèse de récurrence, dim{P^b) — 2{n — 3) — {\V'{a)\ + \V'{b)\) — 
P- 

P^b contient p + 1 points affinement indépendants or Ojj2(7i-3) £ Pâb d'où P^f, 
contient p points linéairement indépendants = {xi^, x^\_q-^^, a;ij, X(Jj_3j 

i e {i,...,p}. 

Soient Xi , Xp les p points de R^^"^^) définis par : 

Xi = {Xia, Xfl^_g^^, ^^0^ , X^j, X^^^gjj,, ^^0^ ) i £ {1, 

Les Xi sont linéairement indépendants par construction et appartiennent à Pab- 
Considérons les partitions : 

Hi = {{a, Vo}, {b},V\ {a, b, vo}} et Ha = {{a}, {b,vo},V\ {a, b, vo}}. 
Ui et 112 vérifient l[î) et l|2ll car uo ^ V(a) U V{b). Soient Xp+i et Xp+2 les 
vecteurs d'incidence associés respectivement à IIi et Ha. Xp+i et Xp+2 ap- 
partiennent à Pab- D'où, {Ojf2(n-2) , Xi, ...,Xp+2} est une famille de p + 3 
points affinement indépendants. Par conséquent, dim{Pab) >p + 2etp + 2 = 
2(n - 2) - {\V'{a)\ + \V'{b)\) = 2{n - 2) - {\V{a)\ + \V{b)\).LsL dernière 
égalité vient du fait que V{a) = V'{a) et V{b) = V'{b). On a donc dim{Pab) ~ 
2{n~2)~{\V{a)\ + \V{b)\). 
. Cas 2 : «0 £ V{a) U V{b), \fvoeV\{a, b}. 

o cas 2.1 : 3vo eV\{a, b} tel que vo ^ V{b). 
De même que dans le cas 1, on peut construire p points Xi linéairement indépen- 
dants appartenant à Pab, ave p = 2(n - 3) - (|1/'(j4)| + Soit Xp+i 
le vecteur d'incidence du ab-séparateur {{a, vo}, {b}, V \ {a, b, vo}}. La famille 
{0]g.2(n-2) , Xi, Xp+i} est une famille de p + 2 points affinement indépendants 
dePai,. On a donc dîm(Pa&) > p+letp+1 = 2(n- 3) - (|F'(a)| + 1 V'(&)|) + 
1 = 2(n- 2) - {\V{a)\ + \V{b)\) car \V{a)\ = \V'{a)\ + 1 et \V{b)\ = \V'{b)\. 
On a donc dim{Pab) = 2(n - 2) - {\V{a)\ + \V{b)\). 
Par symétrie, on obtient un résultat identique si vo ^ A. 

o cas 2.2 : V \ {a, b} — V{a) r]V{b), i.e. tout point de V \ {a, b} est adjacent 
à a et à fe. On a alors dim{Pab) = or j V(a)j = n — 2 et j V(b)| = n — 2 d'où 
2(n- 2) - {\V{a)\ + \V{b)\) = et dim{Pab) = 2{n-2) - {\V{a)\ + 
□ 

2.2 Inégalités valides 

Dans cette partie, nous présentons tout d'abord de nouvelles inégalités valides 
basées sur la structure des chaînes entre a et b. Nous nous intéressons ensuite 
aux inégalités associées aux sous-graphes connexes de G. 
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2.2.1 Inégalités de chaînes 
Proposition 2.2 

Soit Fab une chaîne entre a et h et liVab) l'ensemble des sommets internes à Tab- 
L'inégalité ^ ^ {xia + Xib) < [/(Fab)! — 1 est valide pour Pab- 

A C B 



Preuve : La validité de cette inégalité vient du fait que toute chaîne entre a et 6 
contient au moins un sommet appartenant à C avec {A, B, C} un ab-séparateur 
tel que a e ^ et 6 e _B. 

Pour toute paire (i, j ) de sommets non adjacents de V, soit Qij le nombre max- 
imum de chaînes sommets-disjoints entre i et j. 

On déduit immédiatement de la proposition 12. 21 que pour tout {A,B,C} un ab- 
séparateur tel que a £ A&tb £ B,\e cardinal de C est minoré par aat- 

Pour chaque paire de sommets i et j non adjacents de 1/ \ {a, fe}, les inégalités 
suivantes sont valides pour Pab : 

n-2 

^(Xfca + Xkb) - Ctij * (2 - Xia - Xjb) < U - Uij (11) 

ft = l 

n-2 

^(Xfca + Xkb) - Ctij * (2 - Xja - Xtb) < TL - Oij (12) 
fc = l 

Les inégalités l ll It et l ll2t traduisent le fait que si {A, B,C} est un ab-séparateur 
tel que l'un des sommets i et j est dans A et l'autre dans B, alors \C\ est supérieur 
à Oij. 

2.2.2 Inégalités associées aux sous-graphes connexes de G = {V, E) 

Soit V' un sous-ensemble de V tel que le graphe G' induit par V' est connexe et 
que|l^'| >/3(n). 

Soit = Min{âij, i € V' , j € V' , ^ E} où âij est le nombre 

maximum de chaînes sommets-disjoints entre i et j dans G'. 

^Obtention de ciij: on construit un réseau de transport, avec i et j pour entrée et sortie, pour lequel un 
flot cj) de valeur Oij maximum entre i et j détermine o^j chaînes sommets-disjoints de G = {V, E) entre 
i et J. Pour cela, on remplace dans G chaque arête par deux arcs de dù'ections opposées puis on dédouble 
tout sommet x distinct e i et j par deux sommets x' et x" reliés par un arc (x', x") de capacité 1. Tout 
arc qui allait vers x est remplacé par un arc de capacité infinie allant vers x' et tout arc qui sortait de x 
est remplacé pai' un arc de capacité infinie sortant de x" . Le flot maximum entre i et j représente alors le 
nombre maximum de chaînes sommets-disjoints entre î et j. 

Pour chaque paire de sommets non adjacents, le flot maximum a^j est obtenu à l'aide du code H1_PR 
d'Andrew Goldberg disponible à l'adresse http://www.avglab.com/andrew/soft.html Ce 
code est une version plus robuste du code PRF disponible à la même adresse. Il s'agit d'une implémen- 
tation de la méthode Push-Relabel pour les problèmes de flot maximum / coupe minimum (voir | 8 | |9 |) qui 
est à notre connaissance la meilleure méthode actuellement disponible pour résoudre les problèmes de flot 
maximum. 
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Proposition 2.3 

L'inégalité {xia + Xit) < \V'\ — Min{aQ , \V'\ — est valide pour 

Pab- 

Cette inégalité est intéressante dans le cas de graphes pour lesquels 
(Xmin = Min{aij, i,j £ V, ^ E} est petit mais dont certains sous- 

ensembles de sommets V induisent des sous-graphes connexes pour lesquels «q 
est grand. 



3 Modèle sur les arêtes 
3.1 Formulation 

Le problème est formulé ici en utilisant des variables sur les arêtes du graphe G = 
{V, E). On cherche donc un sous-ensemble d'arêtes F Ç i5 tel que F = FiVJ F2 
avec i^i n F2 = 0, pour lequel il existe une partition {A, B,C} vérifiant l|T) et l|2), 
aG A,bG B. 

F est un sous-ensemble de E arêtes-séparateur. 

Le sous-graphe (Vi U V2,F) induit par F, où Vi = A U S et V2 — C, est un 
graphe biparti. Les résultats classiques pourront donc être utilisés ici. 

Le polyèdre associé au problème des sous-ensembles arêtes-séparateurs est : 
Qab = Conv{x £ R^, 3 un arêtes-séparateur dont X sst le vecteur d'incidence }. 

Pour toute chaîne F , x(r) = ^^x(c)- On note Fat une chaîne entre a et b. 

Si X 6st vecteur d'incidence d'un arêtes-séparateur F, alors x vérifie les inégal- 
ités suivantes : 

x(ra6)>2 vr„6 

(13) 

X(r) -x(ra6\r) < |r| - 1 Vr„6 avecV ÇVat et \r\ est impair 

(14) 

x(r) — x{^ \ r) < |r| — 1 v $ cycle avec F Ç $ et |F| est impair 

(15) 

x(e)G{0,l} VeGÊ 

(16) 

Les inégalités l ll3t et l ll4t résultent de l'appartenance de a et fo à la même classe 
d'un graphe biparti : toutes les chaînes les reliant dans ce graphe sont paires et de 
longueur supérieure ou égale à 2. L'inégalité l ll5t traduit l'absence de cycle impair 
dans un graphe biparti. 






Fl \ / F2 

F = i^i U F2 un arêtes-séparateur 



Proposition 3.1 

Qab = Conv{x G : X vérifie à ^6^} 
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Soit Qab le polyèdre défini par : 

Qab = Conv{x ■■ X vérifie (Ï3}à0}etO< x(e) < 1, Ve G E}. 
Qab est le polyèdre issu de la relaxation du système précédent. 



IMPLÉMENTATION ET RÉSULTATS EXPÉRIMENTAUX 



On a immédiatement Qab Ç Qab mais Qab 7^ Qat- En effet, l'exemple ci- 
dessous exhibe un point extrême fractionnaire de Qat '■ 



Ce point viole l'inégalité 

X(ei) + x(e4) + x(e8) + x(ei2) > x(e2) + xi^s) + xi^a) + xi^?) - 1 
qui est valide pour Qat- 



4 Implémentation et résultats expérimentaux 
4.1 Modèle 



Données : 

• G = (V, -B) un graphe connexe non oriente 

• f3{n) entier tel que 1 < (3{n) < n, (ici, f3{n) = ^ ) 

• a,b affectés respectivement a Aet B (ici a, b sont des sommets fictifs) 

Min{a.,, i(iV,j£V, (ij) i E} 

Sans perte de généralité, on peut supposer que \A\ < \B\. Comme on a également 
1^1 + l-B| < n — Qmm, il en résulte < [-^^^^^I^^J , ce qui est traduit par ([22ll 
dans notre modélisation. 




a 
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Modélisation : 

n 

Maximiser ^(a;^ + Xit) 

i = l 

S.C. : 

x,a G {0, 1}, Vî e V (17) 

x,a + xa>l, Vie y (18) 

x^a + Xjb>l: 'i{i,j)£E (19) 

Xja + xa>l, y{i,j)€E (20) 

n 

(21) 

n 

1<E^"'^ L "~^""" J (22) 

n n 

l<Y^Xii< P{n) , 1 < ^ < I3{n), (23) 



4.2 Résultats 

Nous avons travaillé exclusivement sur les instances étudiées par Balas et De 
Souza pour disposer d'éléments de comparaison. Notre programme linéaire mixte 
a été résolu en utilisant les APIs orientées objet Concert Technology pour C++ du 
solveur Ilog Cplex 9.0. Tous les tests ont été réalisés sur un pc portable équipé 
d'un processeur Intel Pentium M 740 (fréquence 1,73 GHz) et de 1GB de RAM 
(DDR2). Balas et De Souza ont utilisé le solveur XPRESS dont les performances 
sont comparables à celles d'Ilog Cplex. En revanche, leurs résultats sont obtenus 
sur un pc de bureau plus puissant que le nôtre (Intel Pentium 4 à 2,5 GHz et 2GB 
de RAM). 

Tous les graphes ont moins de 150 sommets excepté myciel7 qui en a 191. 
Pour chaque série d'instances, un tableau présenté en annexe rassemble les carac- 
téristiques des graphes qui les composent. Ce tableau fournit le nombre de som- 
mets de chaque graphe, le nombre d'arêtes, la densité, la borne l3{n) associée et la 
valeur de qui donne une première indication sur l'efficacité éventuelle de la 

minoration du cardinal de C par amin. 

4.2.1 Instances DIMACS 

Les instances de cette catégorie sont issues du challenge DIMACS sur la coloration 
de graphe et disponibles à l'adresse : 

fhttp : / / mat . gsia . cmu . edu/COLOR/instances ■ htmTI 
Nos résultats sont présentés dans les tableaux qui suivent. La deuxième colonne in- 
dique le nombre de noeuds que nous avons obtenus, la troisième donnant la valeur 
correspondante pour Balas et De Souza. Les quatrième et cinquième colonnes in- 
diquent respectivement le temps nécessaire à la résolution de l'instance par notre 
algorithme et par celui de Balas et De Souza. L'ajout d'un astérisque précise que 
l'instance n'a pas été résolue à l'optimale par Balas et De Souza dans le temps 
imparti (1800 secondes). Balas et De Souza ont utilisé plusieurs méthodes, nous 
avons systématiquement rapporté leur meilleur résultat, le critère de choix étant le 
temps de calcul. 
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Table 1: Résultats pour les instances DIMACS 



Instance 


nb. noeuds 


nb. noeuds B. S. 


temps(s) 


temps(s) B.S. 


david 


113 


52 


0.984 


0.19 


DSJC125.1 


132000 


218098 


(1800) 


(1800) * 


DSJC125.1 


523328 


218098 


9783.080 


(1800) * 


DSJC125.5 


120 


16768 


8.610 


(1800) 


DSJC125.9 





34241 


0.703 


794.36 


games 1 20 


8088 


96920 


121.047 


429.02 


miles500 


145 


318 


7.156 


2.11 


niiles750 


938 


369 


62.797 


9.83 


miles 1000 


168 


65 


19.062 


13.62 


myciel3 





19 


0.000 


0.00 


niyciel4 


37 


48 


0.109 


0.03 


niyciel5 


290 


169 


0.828 


0.28 


niyciel6 


551 


458 


12.094 


5.14 


mycielV 


14737 


2441 


880.875 


160.59 


queen6_6 





81 


0.016 


1.42 


queen7_7 





263 


0.063 


7.78 


queen8_8 





3533 


0.031 


42.44 


queen8_12 


3175 


162911 


77.297 


(1800) * 


queen9_9 


959 


291471 


14.375 


1067.5 


queenlO_10 


5412 


126778 


105.765 


(1800) * 


queenll_ll 


13410 


55220 


456.781 


(1800) * 


queenl2_12 


24635 


25897 


1245.080 


(1800) * 


moyenne 


9751 


49339 


229.223 


634.97 



En comparant nos résultats à ceux de Balas et De Souza sur les instances DI- 
MACS, nous obtenons : 

• un nombre moyen de noeuds de 9751 contre 49339. 

• un temps moyen de résolution de 229,223 secondes contre 634,97. 

• 20 solutions optimales (en moins de 30 minutes) contre 15. 

On constate en particulier une très bonne efficacité de notre méthode sur les in- 
stances de types n-reines, pour lesquelles amin est élevé, ainsi que sur DS JC 1 2 5 . 5 
(a„i„=5l) et DS JC125 . 9 (ami„=l03) résolues à l'optimal respectivement en 8,61s 
(120 noeuds) et 0,703s (0 noeuds) alors que la première n'est pas résolue en moins de 
30 minutes par Balas et De Souza et que la seconde est résolue en 794,36s. Pour 
l'instance D S JC 1 2 5 . 1 , la meilleure valeur entière obtenue par Balas et De Souza 
en 30 minutes est 89 avec une borne supérieure à 102,16. Nous obtenons dans 
le même temps une meilleure valeur entière de 91 avec une borne supérieure à 
101. Nous obtenons une solution optimale en 9783,080 secondes et la valeur de 
l'objectif est alors de 91, soit notre meilleure valeur entière après 30 minutes de 
calcul. 



4.2.2 Instances MatrixMarket 

Cette catégorie d'instances contient les graphes d'intersection de matrices de sys- 
tèmes linéaires de type Ax = b. Le graphe d'intersection d'une matrice A, noté 
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G (A), a un sommet pour chaque colonne de A et une arête entre une paire de som- 
mets si le produit scalaire des colonnes correspondantes est non nul. L'existence 
d'une arête dans G{A) signifie donc qu'il existe une équation dans Ax — h 
dans laquelle les deux variables Xi et Xj ont des coefficients non nuls. Ainsi, si 
le système linéaire est résolu par une méthode de type diviser pour résoudre {di- 
vide and conquer ), il pourra être divisé en deux sous-systèmes plus petits résolus 
séparément. La solution complète de ces sous-systèmes dépend des valeur des 
variables appartenant aux deux sous-systèmes. Le coût algorithmique de la fusion 
des sous-systèmes nécessaire à l'obtention de la solution du système original aug- 
mente avec le nombre de variables communes aux sous-systèmes. L'efficacité de 
l'algorithme impose aussi que la taille des sous-systèmes soit bornée par une frac- 
tion de la taille du système initial. Le problème du choix de la meilleure façon de 
partitionner le système linéaire se ramène donc à un VSP défini sur G{A). Les 
matrices utilisées pour construire ces instances peuvent être obtenues à l'adresse 
[http : / / math . nist . gov/MatrixMarket Elles correspondent à des prob- 
lèmes pratiques issus de la physique, la météo, l'économie, etc. A partir de ces ma- 
trices, Balas et De Souza ont construit trois catégories d'instances, MM-I, MM-II 
et MM-HD, selon les propriétés des matrices d'origine (voir [2]). Elles peuvent 
être obtenues à l'adresse : 

|http : / / www ■ ic ■ uni camp . br/~cid/Problem- instances /VSP ■ html| 
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Les graphes de la catégorie MM-I ont au moins 20 et au plus 100 sommets. 
Leur densité varie de 0, 06 à 0, 67 



Table 2: Résultats pour les instances MatrixMarket MM_1 



Instance 


iioeiiils 


/;/). iiociuh C.V. 


tciiips(s) 


iciiip.st s) B.S. 


ash219 


137 


111 


1.187 


18 


dwt72 




Jy 


U. JUU 


07 




57 


119 


0.500 


15 


UW 1>JVJ 




31 


0.312 


06 




49 


51 


0.484 


06 


dwt 59 


137 


201 


0.969 


30 




13 


13 


0.125 


02 


ash85 


593 


971 


8.610 


2 87 


HwtR? 


109 


241 


3.031 


90 


llllJJ w \J1 u 


428 


380 


1.532 


59 




195 


185 


1.578 


51 


wilIS? 


55 


17 


0.531 


05 


can96 


7787 


259231 


154.921 


1131.60 


steam3 


108 


313 


1.828 


0.81 


curtis54 


77 


167 


0.625 


0.20 


can73 


1138 


6997 


8.266 


21.82 


bfw62a 


7 


11 


0.579 


0.05 


ibm32 


69 


72 


0.313 


0.07 


porcs_l 


48 




0.204 


0.09 


canôl 


90 


5 


2.109 


0.82 


bcsstkOl 





39 


0.063 


1.63 


can24 





17 


0.000 


0.05 


fidapmOS 


44 


5 


0.422 


0.10 


fidap005 


28 


5 


0.109 


0.04 


moyenne 


467 


11219 


7.867 


48.46 



En comparant nos résultats à ceux de Balas et De Souza sur les instances MM- 
L nous obtenons : 

• un nombre moyen de 467 noeuds contre 11219. 

• un temps moyen de résolution de 7,867 secondes contre 48,46. 
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Les graphes de la catégorie MM-n ont au moins 104 et au plus 147 sonmiets. 
Ils ont été obtenus à partir de matrices dont le nombre de colonnes varie entre 100 
et 200. Le préfixe L12 5 est ajouté au nom de la matrice initiale lorsque le graphe 
est issu de la sous-matrice formée de ses 125 premières colonnes. La densité de 
ces graphes varie de 0, 05 à 0, 93. En comparant nos résultats à ceux de Balas et 



Table 3: Résultats pour les instances MatrixMarket MM_n 



Instance 


nb. noeuds 


nb. noeuds B.S. 


temps(s) 


temps(s) B.S. 


L125.ash608 


142 


82 


2.375 


0.50 


L125.willl99 


135 


119 


2.562 


0.35 


L125.west0167 


69 


73 


1.906 


0.22 


ash331 


183 


212 


2.812 


0.82 


wcslO 1 32 


105 


85 


4.297 


0.42 


rwl36 


1635 


17386 


28.875 


57.57 


bcspwr03 


99 


73 


2.422 


0.38 


gre 1 15 


5573 


6659 


82.812 


28.77 


L125.dw_162 


469 


214 


8.360 


1.15 


L125.can_187 


786 


5353 


16.984 


23.22 


L125.gre_185 


1300 


2337 


30.578 


31.40 


L125.can_161 


15063 


218627 


580.797 


(1800) * 


L125.1opl63 


308 


6319 


16.031 


33.22 


can 144 


1362 


64203 


31.250 


443.49 


lund_a 


385 


2705 


45.891 


155.29 


L125.bcsstk05 


331 


2301 


30.812 


54.29 


L125.dwt_193 


259 


129 


29.719 


73.74 


L125.fs_183_l 


240 


1175 


35.016 


36.21 


bcsstk04 


600 


89 


26.796 


80.01 


arc 130 


268 


129 


99.328 


137.42 


moyenne 


1466 


16414 


53.981 


147.92 



De Souza sur les instances MM-II, nous obtenons : 

• un nombre moyen de 1466 noeuds contre 16414. 

• un temps moyen de résolution de 53,981 secondes contre 147,92. 

• toutes nos solutions sont optimales (en moins de 30 minutes) y compris pour 
L 1 2 5 . can 1 6 1 non résolue par Balas et De Souza. 



Les graphes de la catégorie MM-HD ont tous une densité supérieure à 0, 35. Ils 
contiennent 80, 100 ou 120 sommets. Les préfixes L80, LlOO et L120sont 
ajoutés au nom de la matrice initiale lorsque le graphe est issu de la sous-matrice 
formée de ses 80, 100 ou 120 premières colonnes. 
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Table 4: Résultats pour les instances MatrixMarket MM-HD 



Instance 


nb. noeuds 


nb. noeuds B.S. 


temps(s) 


temps(s) B.S. 


L100.steam2 


175 


167 


8.781 


30.98 


LlOO.cavityOl 


115 


7 


10.657 


3.52 


LSO.cavityOl 


125 


15 


4.578 


2.22 


L80.fidap025 


52 


7 


2.625 


1.14 


L120.fidap025 


82 


21 


9.61 


9.27 


L80.steam2 


115 


29 


4.218 


3.15 


L100.fidap021 


119 


5 


5.469 


2.64 


L100.fidap025 


310 


11 


13.578 


5.41 


L120.cavity01 


138 


15 


16.375 


16.57 


L120.fidap021 


192 


27 


16.765 


17.36 


L80.fidap021 


99 


7 


4.031 


2.43 


L120.rbs480a 


209 


233 


21.282 


98.84 


L120.wm2 


171 


33 


18.781 


19.27 


L100.rbs480a 


156 


63 


13.953 


11.91 


L80.wm3 


184 


13 


8.406 


3.81 


LSO.wml 


155 


59 


8.64 


11.13 


L80.rbs480a 


43 


5 


3.5 


0.97 


L80.wm2 


153 


15 


7.36 


2.89 


LlOO.wmB 


179 


9 


19.656 


4.79 


L120.e05r0000 


82 


11 


9.828 


7.51 


LlOO.wml 


196 


27 


27.782 


11.65 


L120.fidap022 


219 


177 


30.782 


94.06 


L100.wm2 


190 


15 


18.485 


7.37 


L100.fidapm02 


255 


13 


25.125 


6.35 


L120.fidap001 


237 


27 


37.532 


22.56 


L100.e05rOOOO 


120 


47 


7.937 


14.94 


L120.fidapm02 


223 


11 


42.953 


12.29 


L80.fidapm02 


50 


7 


3.281 


1.97 


LlOO.fidapOOl 


100 


39 


5.141 


9.82 


Ll()0.1idap()22 


168 


143 


13.985 


30.62 


L80.c05rOOOO 


133 


3 


7.937 


0.96 


L80.fidap001 


84 


1 


3.14 


0.98 


L80.fidap022 


75 


171 


2.172 


10.05 


L80.fidap002 


83 


9 


3.328 


1.71 


L80.fidap027 


81 


3 


4.469 


1.18 


L100.fidap027 


187 


5 


23.437 


4.15 


L100.fidap002 


102 


3 


6.234 


2.38 


L120.fldap002 


121 


71 


10.453 


37.86 


L120.fidap027 


229 


9 


50.719 


9.07 


moyenne 


146 


39 


13.666 


13.74 



En comparant nos résultats à ceux de Balas et De Souza sur les instances 
MM-HD, nous obtenons un nombre moyen de 146 noeuds contre 39 et un temps 
moyen de résolution de 13,666 secondes contre 13,74. 
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4.3 Synthèse 

Il semble à l'issue de ces tests que notre algorithme soit plus performant que celui 
de Balas et De Souza sur les instances difficiles. Ses performances sont directe- 
ment liées à la structure du graphe initial. En effet, la séparation d'un graphe dans 
lequel toute paire de sommets non adjacents est associée à un nombre élevé de 
chaînes sommets-disjoints est extrêmement rapide pour les instances testées. Les 
exemples de DSJC125.5, DSJC125.9, queenlO_10, queen9_9 mettent en évidence 
ce constat. Sur ces quatre instances pour lesquelles amin est élevé, nous explorons 
en moyenne 1623 noeuds contre 117314 pour Balas et De Souza. Notre temps de 
résolution moyen est de 32,363 secondes et toutes nos solutions sont optimales. 
Balas et De Souza résolvent seulement deux de ces instances à l'optimal en 930,93 
secondes en moyenne. La densité du graphe initial affecte les performances de 
l'algorithme de Balas et De Souza. Notre algorithme est pour sa part mieux adapté 
au traitement des graphes denses, susceptibles d'être dotés de valeurs élevées pour 

5 Conclusions et perspectives 

Dans ce rapport, nous étudions une formulation du problème du séparateur d'un 
graphe G = {V,E) connexe non orienté grâce à un programme linéaire mixte, 
donné par Balas et De Souza, que nous avons enrichi. Nos travaux nous ont per- 
mis d'obtenir la dimension du polyèdre Pat associé aux ab-séparateurs du graphe, 
enveloppe convexe de l'ensemble des vecteurs d'incidence des partitions de V véri- 
fiant le VSP avec a et b fixés. Nous avons introduit de nouvelles inégalités valides 
en considérant les chaînes sommets-disjoints entre a et b. L'algorithme que nous 
avons implémenté utilise ces inégalités et nous permet d'améliorer sensiblement 
les résultats obtenus par Balas et De Souza |2|. 

Nous proposons parallèlement un modèle introduisant des variables sur les arêtes 
du graphe. 

Nous poursuivons actuellement notre étude du VSP dans une direction théorique 
et dans une direction expérimentale. Du point de vue théorique, l'approche polyé- 
drale du problème nous conduit à rechercher les facettes définies par les inégalités 
valides de notre modélisation. Nous souhaitons également compléter le modèle 
utilisant des variables sur les sommets par l'ajout de celui utilisant des variables 
sur les arêtes pour poursuivre la recherche de nouvelles inégalités valides amélio- 
rant la description polyédrale du problème. Du point de vue expérimental, nous 
travaillons à la mise au point de stratégies de séparation des inégalités que nous 
avons introduites. Nous préparons également une série de tests qui nous permettra 
d'évaluer l'efficacité de notre algorithme sur des instances de grandes tailles. 
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Annexes : 

Ces annexes rassemblent les caractéristiques des graphes qui composent les séries 
d'instances sur lesquelles nous avons travaillé. Dans chaque tableau, n est le nom- 
bre de sommets, e est le nombre d'arêtes, d est la densité du graphe (d = 
/3(n) est la borne sur le cardinal des ensembles AtiB il3{n) = ^), amin est le 
nombre minimum de sommets appartenant au séparateur C. 

A Instances DIMACS 



Table 5: Instances DIMACS 



Instance 


n 


e 


a 


p[n) 




t^min / ^ 


david 


87 


406 


0.11 


58 


1 


0.011 


DSJC 125.1 


125 


736 


0.09 


83 


5 


0.040 


DSJC125.5 


125 


3891 


0.50 


83 


51 


0.408 


DSJC125.9 


125 


6961 


0.90 


83 


103 


0.824 


gamesl20 


120 


638 


0.09 


80 


2 


0.017 


milesSOO 


128 


1170 


0.14 


85 


2 


0.016 


miles750 


128 


2213 


0.26 


85 


6 


0.047 


mileslOOO 


128 


3216 


0.40 


85 


11 


0.086 


mycieO 


11 


20 


0.36 


7 


3 


0.273 


myciel4 


23 


71 


0.28 


15 


4 


0.174 


myciel5 


47 


236 


0.22 


31 


5 


0.106 


mycielô 


95 


755 


0.17 


63 


6 


0.063 


mycielV 


191 


2360 


0.13 


127 


7 


0.037 


queen6_ 6 


36 


290 


0.46 


24 


15 


0.417 


quccn7_ 7 


49 


476 


0.40 


32 


18 


0.367 


quccn8_ 8 


64 


728 


0.36 


42 


21 


0.328 


quccn8_ 12 


96 


1368 


0.30 


64 


25 


0.260 


queen9_ 9 


81 


1056 


0.33 


54 


24 


0.296 


queenlO_ 10 


100 


1470 


0.30 


66 


27 


0.270 


queenll_ 11 


121 


1980 


0.27 


80 


30 


0.248 


queenl2_ 12 


144 


2596 


0.25 


96 


33 


0.229 
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B Instances MatrixMarket 



Table 6: Instances MatrixMarket MM-I 



Instance 


n 


e 


d 


P{n) 




amin/n 


ash219 


85 


219 


0.06 


56 


2 


0.024 


dwt72 


72 


170 


0.07 


48 


2 


0.028 


can62 


62 


210 


0.11 


41 


2 


0.032 


dwt66 


66 


255 


0.12 


44 


4 


0.061 


bcspwr02 


49 


177 


0.15 


32 


2 


0.041 


dwt_59 


59 


256 


0.15 


39 


3 


0.051 


bcspwrOl 


39 


118 


0.16 


26 


2 


0.051 


ash85 


85 


616 


0.17 


56 


6 


0.071 


dwt87 


87 


726 


0.19 


58 


4 


0.046 


impcol_b 


59 


329 


0.19 


39 


3 


0.051 


west0067 


67 


411 


0.19 


44 


3 


0.045 


will57 


57 


304 


0.19 


38 


2 


0.035 


can96 


96 


912 


0.20 


64 


14 


0.146 


steam3 


80 


712 


0.23 


53 


8 


0.100 


curtis54 


54 


337 


0.24 


36 


5 


0.093 


can73 


73 


652 


0.25 


48 


10 


0.137 


bfw62a 


62 


639 


0.34 


41 


3 


0.048 


ibm32 


32 


179 


0.36 


21 


4 


0.125 


pores_l 


30 


179 


0.41 


20 


5 


0.167 


canôl 


61 


866 


0.47 


40 


12 


0.197 


bcsstkOl 


48 


622 


0.55 


32 


18 


0.375 


can24 


24 


156 


0.57 


16 


8 


0.333 


fidapm05 


42 


528 


0.61 


28 


12 


0.286 


fidapOOS 


27 


234 


0.67 


18 


9 


0.333 
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Table 7: Instances MatrixMarket MM_II 



Instance 


n 


e 


d 


Pin) 






L125.ash6()8 


125 


390 


0.05 


83 




0.016 


L125.willl99 


125 


386 


0.05 


83 


1 


0.008 


L125.west0167 


125 


444 


0.06 


83 


1 


0.008 


ash331 


104 


331 


0.06 


69 


2 


0.019 


west0132 


132 


560 


0.06 


88 


1 


0.008 


rwl36 


136 


641 


0.07 


90 


1 


0.007 


bcspwr03 


118 


576 


0.08 


78 


2 


0.017 


gre_115 


115 


576 


0.09 


76 


3 


0.026 


L125.dw_162 


125 


943 


0.12 


83 


6 


0.048 


L125.can_187 


125 


1022 


0.13 


83 


12 


0.096 


L125.gre_185 


125 


1177 


0.15 


83 


6 


0.048 


L125.can_161 


125 


1257 


0.16 


83 


12 


0.096 


L125.1opl63 


125 


1218 


0.16 


83 


9 


0.072 


can 144 


144 


1656 


0.16 


96 


18 


0.125 


lund_a 


147 


2837 


0.26 


98 


11 


0.075 


L125.bcsstk05 


125 


2701 


0.35 


83 


17 


0.136 


L125.dwt_193 


125 


2982 


0.38 


83 


21 


0.168 


L125.fs_183_l 


125 


3392 


0.44 


83 


1 


0.008 


bcsstk04 


132 


5918 


0.68 


88 


48 


0.364 


arcl30 


130 


7763 


0.93 


86 


21 


0.162 
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Table 8: Instances MatrixMarket MM-HD 



Instance 


71 


c 


d 


0in) 






ijiUU.SlcamZ 


1 nn 

lUU 


1 /oo 


n '^fi 

U.JO 


AA 


1 A 
10 


n 1 An 
u.iou 


L, i UU.CaVliyU 1 


1 nn 

lUU 


1 OAzl 


n '^7 


AA 
00 



o 


n non 

U.UoU 


L,ou.caviiyui 


on 
Î5U 


1 oni 
IZUi 


n 'îo 
U. Jo 


^T 



O 


n 1 nn 
U.lUU 


J_,oU.nQapUZj 


on 

oU 


1 oni 

IZUi 


n '^0 


^T 


1 

IZ 


n 1 ^n 

U. 13U 


ijizu.naapuzD 


1 on 
izu 


0707 
Z / o / 


n '^Q 
u. 


on 

oU 


1 
lo 


n 1 ^n 
u.i ou 


LoU.sieamz 


on 


1 070 

IZ /Z 


n /in 

U.4U 


^T 


1 
IZ 


n 1 <n 

U.IDU 


T 1 nn firlnr^n') 1 


1 nn 

lUU 


onoo 

ZUZo 


n zi 1 


AA 
00 




n 1 '>n 

U. 1 jU 


T 1 nn fi/^^r>nO'^ 


1 nn 

1 WU 


om 1 

ZUJ 1 


n /i 1 

U.4 1 


AA 
00 


1 

1 z 


n 1 on 

u. 1 ZV/ 


T I m .-..1, ; 1 1 n 1 
L 1 ZU.CcU 11) U 1 


1 on 
1 ZU 


"»070 

zv / z 


n 1 

U.4Z 


S'A 


Q 

o 


u.uo / 


L, 1 zu. Tiuapuz 1 


1 on 
izu 


^n^o 


n /n 


on 
oU 


1 c 

IJ 


n 1 0"^ 
u. 1 ZD 


T sn firiQi-*no 1 
j_,ou.naapuzi 


on 

oU 


1 ^Al 
IJ^ 1 


n 




1 1 
1 1 


n 1 To 
u. 1 


Lizu.rDS4oUa 


1 on 
izu 




n /lA 
U.40 


on 
oU 


OT 

Zj 


n 1 oo 
u.iyz 


T 1 Of\ \-i7rv\0 

ijizu.wniz 


1 on 

izu 




n /i7 
U.4 / 


on 


1 
1 


n nno 
U.UUo 


T 1 nn t-v»o4sna 

J_<iUU.rDS4oUa 


1 nn 

lUU 




n ^0 

U. jZ 


AA 
00 


OA 

zo 


n OAn 
u.zou 




on 
oU 


1 7TQ 


n 

U. j J 


^T 


1 
1 


n ni T 
U.Ul j 


T on ïirmi 


on 
bU 


1 70A 
1 /oO 


n ^7 
U. j / 


^T 


z 


n no^ 
U.UZj 


T sn t-KcAiîno 


on 

oU 


1 1 Q 


n 

U. Jo 


^T 


1 
lo 


n 00^ 

U.ZZJ 


T on \\7rrtl 


on 

oU 


1 OAO 
lo^o 


n 

U. Jo 


^T 


1 
1 


n ni T 

U.Ul j 


T 1 nn ii7mQ 


1 nn 

lUU 




n 


AA 
00 


1 
1 


n m n 
U.UlU 


T 1 on »an^rnnnn 

L i ZU . eu jrUUUU 


1 on 
izu 


Al 77 


n 
u. 


on 

oU 


^n 

:>u 


n o^n 

U.ZjU 


T 1 nn \\?m 1 
ijiuu.wnii 


1 nn 

lUU 


OQ^A 
ZVjO 


u.ou 


AA 
00 


1 


n ni n 
u.uiu 


T 1 on ■fi/-ioi-»noo 
zu. naapuzz 


1 on 
izu 


A ^n7 

4 jU / 


n An 
U.dU 


on 
oU 


^n 
:>U 


n o'^n 
U.ZjU 


T 1 nn \ï/mO 
Liuu.wniz 


1 nn 

lUU 




n 1 

U.Ol 


AA 
00 


1 
1 


n ni n 
u.uiu 


T 1 nn firititAmno 
j-< i uu . nuapmuz 


1 nn 

lUU 


^non 
juyu 


n ^^0 
u.oz 


AA 
00 


07 
Z / 


n 07n 
u.z /u 


T 1 on fi^iit^nni 
L i zu. nuapuu 1 


1 on 
izu 


00 


n AT 


on 
oU 




n ooo 
u.zvz 


T 1 00 pOSrOOOO 


1 00 

lUU 




f\A 


UU 


^0 


^00 


L120.fidapm02 


120 


4626 


0.65 


80 


27 


0.225 


L80.lidapm()2 


80 


2058 


0.65 


53 


27 


0.338 


LlOO.fidapOOl 


100 


3379 


0.68 


66 


36 


0.360 


L100.fidap022 


100 


3359 


0.68 


66 


38 


0.380 


L80.e05rOOOO 


80 


2150 


0.68 


53 


19 


0.238 


LSO.fidapOOl 


80 


2288 


0.72 


53 


26 


0.325 


L80.fldap022 


80 


2403 


0.76 


53 


39 


0.488 


L80.fidap002 


80 


2422 


0.77 


53 


27 


0.338 


L80.fidap027 


80 


2518 


0.80 


53 


24 


0.300 


L100.fidap027 


100 


4014 


0.81 


66 


27 


0.270 


L100.fidap002 


100 


4057 


0.82 


66 


34 


0.340 


L120.fidap002 


120 


5886 


0.82 


80 


52 


0.433 


L120.fidap027 


120 


6051 


0.85 


80 


35 


0.292 



